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It is well-known that even high quality cosmic microwave background (CMB) observations are 
not sufficient on their own to determine the equation of state of the dark energy, due to the effect 
of the so-called geometric degeneracy at large multipoles and the cosmic variance at small ones. In 
contrast, we find that CMB data can put tight constraints on another fundamental property of the 
dark energy, namely its coupling to dark matter. We compare the current high-resolution CMB data 
to models of dark energy characterized by an inverse power law or exponential potential and by the 
coupling to dark matter. We determine the curve of degeneracy between the dark energy equation 
of state and the dimensionless Hubble parameter h and show that even an independent perfect 
determination of h may be insufficient to distinguish dark energy from a pure cosmological constant 
with the current dataset. On the other hand, we find that the interaction with dark matter is 
firmly bounded, regardless of the potential. In terms of the dimensionless ratio f3 of the dark energy 
interaction to gravity, we find /3 < 0.16 (95% c.l.). This implies that the effective equation of state 
between equivalence and tracking has been close to the pure matter equation of state within 1% 
and that scalar gravity is at least 40 times weaker than tensor gravity. Further, we show that an 
experiment limited by cosmic variance only, like the Planck mission, can put an upper bound (3 < 
0.05 (95% c.l.). This shows that CMB observations have a strong potentiality not only as a test of 
cosmic kinematics but also as a gravitational probe. 

I. INTRODUCTION 

Four years after the first observational hints about the existence of a dominant component of unclustered matter 
with negative pressure (Riess et al. 1998, Perlmutter et al. 1999), the so-called dark energy or quintessence (Wetterich 
1988; Ratra & Peebles 1988; Frieman et al. 1995; Caldwell et al. 1998), there are still very few indications as to its 
nature. The main reason, perhaps, is that we lack any specific theoretical suggestion on the properties of the dark 
energy, i.e. on its self-interaction potential and on how it interacts with the other cosmological components. At this 
stage, all we can do is to explore a wide range of phenomenological models specified by the potential and the coupling 
to the other fields in order to provide an overall best fit to the current data. 

Several works have tried to constrain dark energy models with the recent high resolution CMB data (Netterfield 
et al. 2001, Lee et al. 2001, Halverson et al. 2001) leading in some cases to interesting but conflicting results. In 
Amendola (2001) the coupling dimensionless parameter (3 that measures the strength of the interaction of dark energy 
to dark matter with respect to the gravitational interaction was constrained to be smaller than 0.1 roughly, adopting 
an exponential potential. The slope of the potential, on the other hand, was found to be essentially unconstrained by 
the data. Baccigalupi et al. (2002) found that, among power-law potentials V ~ 4>~ a and fixing the Hubble constant 
to ft- = 0.65 (in units of 100 km/sec/Mpc), values of a around unity are favored while the case of pure cosmological 
constant, a — 0, was less likely by a factor of five roughly. Corasaniti & Copeland (2001), on the other hand, have 
shown that the pure cosmological constant gives the best fit to the data, particularly for as concerns the position of 
the acoustic peaks. Bean & Melchiorri (2002) also conclude that the pure cosmological constant gives the best fit to 
the CMB data when the prior on the Hubble constant is broadened to h = 0.72 ± 0.08 . All the papers above, and 
the present one, assume a flat spatial curvature. 

In this paper we extend the previous studies in two respects. First, we aim at constraining not only the dark energy 
self-interaction (i.e. its potential) but also its interaction with matter. To our knowledge, this is the first work that 
simultaneously constrain the dark energy potential and coupling. Second, we put only very weak restrictions on the 
cosmological parameters h, Qb,c (the density parameters of baryons and cold dark matter, respectively), and n s (the 
slope of the scalar perturbations). The conflicting results above are in fact due mostly to different priors. 

These generalizations allow us to formulate the main conclusion of the paper: we find that the present CMB data 
are capable to put a strong constraint on the coupling but not on the scalar field potential or equation of state. 
In fact, the degeneracy between h and Wj,, the dark energy equation of state (Huey et al. 1999) almost erases the 
sensitivity of the CMB to the dark energy potential. In sharp contrast, as it will be shown below, the CMB spectra 
are very sensitive to the dark energy coupling, since the latter determines the equation of state for a long stage after 
equivalence and there are no strong degeneracies with the other cosmological parameters. 
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II. DARK-DARK COUPLING 



Let us consider two components, a scalar field <f> and CDM, described by the energy-momentum tensors T^ v ^ and 
^)w(c) t respectively. General covariance requires the conservation of their sum, so that it is possible to consider an 
interaction between dark energy and dark matter such that 



-CT, 



(c)<P;v 



(1) 

(2) 



Such a coupling can be obtained from the conformal transformation of a Brans-Dicke gravity (see e.g. Amendola 
1999) and it has been considered several times in literature starting from Wetterich (1988,1995) and Wands et al. 
(1993). It has been discussed in the context of dark energy models in Amendola (1999,2000) and in Wands & Holden 
(2000), Chimento et al. (2000), Billyard & Coley (2000), Chiba (2001), Albrecht et al. (2001), Esposito-Farese & 
Polarsky (2001). In its conformally related Brans-Dicke form has been studied by Uzan (1999), Chiba (1999), Chen 
& Kamionkowsky (1999), Baccigalupi et al. (2000), Sen & Sen (2001). Theoretical motivations in superstring models 
and in brane cosmology have been proposed recently in Gasperini, Piazza & Veneziano (2002) and Pietroni (2002). 

In the flat conformal FRW metric ds 2 = a 2 {—dT 2 +8 i jdx l dx^) the scalar field and dark matter conservation equations 
are 



2H4> + a 2 U^ = Cp c a 2 , 
p c + 3Hp c = -Cp c 4> 



(3) 
(4) 



(dots refer to conformal time) where H = a/a. Suppose now that the potential U(4>) is specified by the following 
relation 



U' = BU 



N 



(5) 



(the prime refers to derivation with respect to (f) which includes power laws, exponential potentials and a pure 
cosmological constant. In the case of power law potential U = A^T a we have 



N = (1 + a) /a 



(6) 



and B = — aA^ x ^ a \ while for the exponential potential N = 1. We consider only the range N > 1 since for negative 
a there are no asymptotically accelerating models; for N — > oo we recover the pure cosmological constant. Assuming 
that the baryons are not directly coupled to the dark energy (otherwise local gravity experiment would reveal a fifth 
force, see Damour et al. 1990) and that the radiation as well is uncoupled (as it occurs if the coupling is derived 
by a Brans-Dicke Lagrangian, see e.g. Amendola 1999), the system of one Einstein equation and four conservation 
equations (for radiation, 7, baryons, b, CDM , c, and scalar field) can be conveniently written introducing the following 
five variables that generalize Copeland et al. (1997): 



k (p 




(7) 



where k 2 = 8nG (notice that w = dloga/dt i.e. the usual Hubble constant). Notice that x 2 ,y 2 ,z 2 , etc. correspond 
to the density parameter of each component. In terms of the independent variable logo we have then the system: 
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where 

The dimensionless constant /? 2 can be seen as the ratio of the dark energy-dark matter interaction with respect to 
gravity. It can be shown in fact (Damour & Nordtvedt 1993, Wetterich 1995) that the force acting between dark 
matter particles can be described in the Newtonian limit as a renormalized Newton's constant G — G(l + 4/3 2 /3). 
The effect of this interaction on structure growth is discussed in Amendola & Tocchini-Valentini (2002). 



III. CRITICAL POINTS AND TRACKING SOLUTIONS 

The system (||) includes several qualitatively different behaviors, already discussed in Amendola (2000), Tocchini- 
Valentini & Amendola (2002). However, for the range of values that are of cosmological interest, the system passes 
through three distinct phases after equivalence. 

</>MDE. Immediately after equivalence, the system enters a matter dominated epoch with a non-negligible <j> con- 
tribution, that we denote 0MDE, in which the dark energy potential density is negligible while the kinetic energy 
density parameter VLk <t> — x 2 of the scalar field gives a constant contribution to the total density. As will be shown in 
the following, the existence of such an epoch is crucial for the constraints that we will be able to put on the coupling. 
Putting y = and neglecting radiation and baryons we obtain the simplified system 

x' = -^x(l-x 2 ) +(3(1 -x 2 ), 

V = o. (10) 



The point y = 0, x = 2/3/3 is a saddle point solution of the system ( |l0| ) for any N and for \/3\ < v / 3/2. The system 
stays on the 0MDE solution until y starts growing. Along this solution the scale factor expands slower than a pure 
MDE, i.e. as 

a ^ t 4/(6+4^)_ (n) 

Since y = on the ^MDE, its existence is independent of the potential, although it has to verified for each potential 
whether it is a saddle. 

Tracking trajectories. Let us now neglect baryons and radiation and put j3 = 0. The tracking solutions found in 
Steinhardt et al. (1999) assume y/x = p and y 2N - 1 w 2N - 2 = q where p, q are two motion integrals. In the limit 
y 2 , x 2 -C 1 it is easy to show that p' = q J = if 

p 2 = AN -3 (12) 
3p 

q = -2M2N-1) (13) 

The same tracking behavior remains a good approximation for small /3. For N = 1 the tracking solution becomes 
actually a global attractor of the dynamical system, see below. In the other cases, the tracking interpolates between 
the 0MDE and the global attractor. 

Global attractors. In the system (||) for N ^ 1 there exists only the attractor x — 0, y — 1 on which the dark energy 
completely accounts for the matter content. For N = 1 the phase space is much richer, and there are several possible 
global attractors, only two of which accelerated (Amendola 2000). One, for ^ > (— j3 + y/18 + /3 2 )/2 , presents a 
constant non-zero fl c : this attractor actually coincides with the tracking solutions and in fact realizes the condition 
p = 1 and y = —3/(2//) as requested by Eq. ( p"2| ) for (3 = 0. This "stationary attractor" can be accelerated if (3 > 2/x 
and could solve the coincidence problem since fl c oc O^; it has been discussed in detail in Amendola & Tocchini- 
Valentini (2001, 2002). The other accelerated attractor occurs when /J, < (—/3 + v / l8 + P 2 )/2 : in this case there are 
no tracking solutions and the global attractor reduces to x = 0, y = 1 as for N ^ 1. These solutions have already been 
compared to CMB for TV = 1 in Amendola (2001). The inclusion of the baryons modifies the considerations above 
but, as long as they are much smaller than the other components, the qualitative behavior of the system remains the 
same. It is to be noticed that the final attractor, on which the dark energy dominates completely the cosmic fluid, is 
yet to be reached, and therefore the existence of an accelerated epoch at the present depends mostly on the tracking. 
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The existence of the 0MDE saddle and of the tracking solutions is crucial for our analysis. In fact, these two 
epochs guarantee that the equation of state of the scalar field is piece-wise constant through essentially all the post- 
equivalence epoch. In the </MVlDE phase the effective parameter of state w e — Ptot/ptot + 1 and the field equation of 
state w<p = p<p/ p<p + 1 are 

w e = l + ^/3 2 , w -2 (14) 

while during the tracking phase 

w e *l, W ^-1^ = -A— (15) 

(the last relation is approximated and it is actually only an upper bound to the present w<f>; it is more precise if is 
identified with the average equation of state after 0MDE rather than the present equation of state). Therefore, the 
cosmic evolution depends on (3 alone during the </>MDE, and on N alone during the tracking. Since the position of 
the acoustic peaks is related to the equation of state through the angular diameter distance, it appears that the CMB 
is able to put direct constraints on both (3 and N. In Fig. 1 we show a typical trajectory that presents in sequence 
the three epochs discussed above. 

There is however a crucial difference between the </>MDE and a tracking for what concerns here: while the present 
dark energy equation of state, set by the tracking, is degenerated with h for as concerns the CMB spectrum (see 
e.g. Huey et al. 1999; Bean & Melchiorri 2002), the equation of state during the <^MDE is not. In fact, the angular 
diameter distance to the last scattering surface cIa is degenerate along lines h(w<f,) for which 

d A ~ f da[w c a+ (h 2 - uj c )a 4 - 3w o]~ 1/2 = const (16) 

Ja dec 

where lo c = fl c h 2 (although this is exact only for (3 = it remains a good approximation even for small non-zero 
values). Assuming a strong prior on ft- a peak emerges in the likelihood for N but then the result is clearly prior- 
dependent. The same holds true for models in which the equation of state is slowly varying (see e.g. Huey et al. 
1999, Doran et al. 2001). On the other hand, the fact that fl^ ^ at decoupling in coupled models implies that the 
effects of the coupling on the CMB are not due solely to the angular diameter distance, and therefore the geometric 
degeneracy can be broken. This is shown in Fig. 2 in which Ci spectra for various values of j3 (all other parameters 
being equal) are shown: the spectra change both in amplitude and in peak's position. 



IV. CMB CONSTRAINTS 



Our theoretical model depends on three scalar field parameters and four cosmological parameters: 

(3,u,N,n s ,h,uj b ,uj c (17) 

where uj b = flbh 2 and to c = fl c h 2 . We calculate the Ci spectra by a modified CMBFAST (Seljak & Zaldarriaga 1996) 
code that includes the full set of coupled perturbation equations (see Amendola 2000 and Amendola & Tocchini- 
Valentini 2001). The other parameters are set as follows: T cmb = 2.72QK, Y He = 0.24, N v = 3.04, r c = 0. In the 
analysis of Netterfield et al. (2002) r c , the optical depth to Thomson scattering, was also included in the general 
likelihood but, in the flat case, was found to be compatible with zero and to have only a minor effect on the other 
parameters. We did not include the cross-correlation between bandpowers and pointing and beam corrections because 
they are not available for all the experiments; we adopted a pure log-normal likelihood. We used as dataset the COBE 
data analyzed in Bond et al. (2000), and the high resolution data of Boomerang (Netterfield et al. 2002), Maxima 
(Lee et al. 2002), and DASI (Halverson et al. 2002). The overall amplitude and the calibration errors of Boomerang 
(10%) and of Maxima and DASI (4%) have been integrated out analytically. The theoretical spectra have been binned 
as the experimental ones. 

In Amendola (2001) we have shown that the dynamics of the system is insensitive to the sign of (3, since </>MDE, 
tracking and the final accelerated epoch do not depend on it. We will consider only (3 > 0. 

The tracking trajectories are selected for each N, (3 by fixing the initial conditions on x, y, and then varying n until 
the present values of h,uj b ,uj c are found. Therefore, the parameter fi is not a free parameter, but is fixed in function 
of the others. In fact, if there is tracking, all trajectories with equal N,/3, h,uj bl uj c always had the same equation of 
state and matter densities, and therefore generate a similar CMB spectrum, so that it is sufficient to select one of 








Figure 1: Numerical solutions of the system ^ for N = 2, /3 = 0.1, uj c = 0.1, uib = 0.02, h = 0.65 plotted against the redshift. 
Upper panel. Long dashed line: radiation; short dashed line: CDM; unbroken thin line: baryons; unbroken thick line: scalar 
field; dotted line: scalar field potential energy. The horizontal thin line marks the kinetic energy density of 0MDE, reached 
just after equivalence. Lower panel. Thick line: effective equation of state; thin line: dark energy equation of state. The labels 
mark the 0MDE, the tracking (T) and the final attractor (A). 



them. There are in general initial conditions that do not lead to tracking (but still undergo </>MDE). These will be 
discussed in another work. 

In order to compare with the previous analyses we assume uniform priors with the parameters confined in the range 
(3 G(0,0.3), N G (1, 8.5), n s G (0.7,1.3), h£ (0.45,0.9), uj b e (0.005, 0.05), oj c G (0.01, 0.3) . The same age 
constraints (> 10 Gyr) used in most previous analyses is adopted here. A grid of ~ 10,000 multipole CMB spectra 
Ce is used as a database over which we interpolate to produce the likelihood function. Notice that we cannot use 
here the rescaling of the high-£ spectra adopted in Bean & Melchiorri (2002) to speed up calculations, due to the 
non-vanishing scalar field energy density during the 0MDE. 



V. LIKELIHOOD RESULTS 



In Fig. 3 we show the likelihood curves for TV, h marginalizing over the remaining parameters. The contours of the 
likelihood plot follow the expected degeneracy of the angular diameter distance. The likelihood is almost degenerated 
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Figure 2: CMB spectra for increasing values of the coupling constant (the other parameters are N = 1.5, h = 0.65, ujb = 
0.01, uj c = 0.1). Notice that the peaks not only move to the right but also change in height. 



along the dotted curve 

h = 0.94 - 0.56uty + 0.03tt$ (18) 

In the plot we also show the expected degeneracy curve from Eq. ( |l^ ) assuming uj c = 0.1. The residual deviation from 
the expected degeneracy line is due to the age prior (>10 Gyr) that favors small h values. Notice that for h > 0.75 
no upper limit to N can be given with the current CMB data no matter how precise h is, and that only assuming 
h < 0.65 it becomes possible to exclude at 95% c.l. the pure cosmological constant. In Fig. 4 we show the likelihood 
for all parameters, marginalizing in turn over the others. Notice that the likelihood for N flattens for TV > 1.2. We 
find the following constraints at 95% c.l.: 

N > 1.5, /3<0.16 (19) 

(the limit on N corresponds to a < 2). The limit on N is however only a formal one: the likelihood never vanishes 
in the definition domain and even N = 1 (the exponential potential) is only a factor of five less likely than the peak 
and certainly cannot be excluded on this basis. Moreover, the lower bound on N is prior-dependent: allowing smaller 
values of h it would weaken. Clearly, if we adopt narrow priors on h, we can indeed obtain more stringent bounds on 
N, as shown in the same Fig. 4. In contrast, the constraint on (3 does not depend sensibly on the prior on h and the 
likelihood for does vanish at large (3. In place of TV and (3 we can use as well the equation of state during tracking 
and during 0MDE, respectively, as likelihood variables, using Eqs. (14) and (|l5|). Then we obtain at the 95% c.l. 



W ^tracking) < 0.8, 1 < W e[(j)MDE) < 1.01 (20) 



This shows that the effective equation of state during <f)MDE, i.e. between equivalence and tracking, is close to unity 
(as in a pure matter dominated epoch) to within one per cent. The striking difference between the level of the two 



7 



N 

5 2 1.5 1.1 




0.2 0.4 0.6 0.8 

w 



Figure 3: Likelihood contour plots in the space N,h marginalizing over the other parameters at the 68,90 and 95% c.l.. The 
white thick lines refer to the coupled case, the black curves to the uncoupled case. The dotted line is the likelihood degeneracy 
curve, the dashed line is the expected degeneracy curve. Notice that only fixing h smaller than 0.65 it would be possible to 
exclude the cosmological constant at 95% c.l.. 

constraints in (|2p| ) well illustrates the main point of this paper: the CMB is much more sensitive to the dark energy 
coupling than to its potential. 

The other parameters are (we give here for simplicity the mean and the one sigma error, while the limit is at the 
95% c.l.) 

n s = 0.99 ±0.05, u b = 0.023 ± 0.004, oj c = 0.092 ± 0.02, h > 0.62 (21) 

The total dark energy density turns out to be > 0.60 (95% c.l.). It appears that the limits on the cosmological 
parameters n s Wb, u c are almost independent of (3, while a non-zero (3 favors higher h (Fig. 4, panel d). It is interesting 
to compare with the current constraints from the Hubble diagram of the supernovae la, where one obtains quite a 
stronger bound on the equation of state, w<f, < 0.4 or N > 1.75 at the same c.l.. 

Finally, in Fig. 4 (panel b) we plot the likelihood for (3 that an experiment with no calibration uncertainty and 
limited only by cosmic variance, like the Planck mission (de Zotti et al. 1999), can achieve. We find (3 < 0.05 (95% 
c.l.) (using the specifications of the MAP satellite we find (3 < 0.1). 

VI. CONCLUSIONS 

The main result of this paper is that CMB observations are a powerful gravitational probe able to constrain 
scalar gravity. In fact, a dark matter-dark energy interaction would obviously escape any local gravity experiment: 
cosmological observations like the CMB are then the only way to observe such a phenomenon. Since observations 
require the baryons to be decoupled from dark energy (or coupled much more weakly than dark matter), the search 
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Figure 4: Marginalized likelihood for tracking trajectories. The equation of state in panel a is w^, = 1/(2N — 1). In panels a 
and b we plot as a dotted line the likelihood assuming h = 0.65 ± 0.05 and as dashed line h = 0.75 ± 0.05 (gaussian prior). In 
panel b the long-dashed line is the Planck-like likelihood. In panel d the long-dashed line is for j3 — 0. 



for a non-zero (3 is in fact also a test of the equivalence principle. We found that current CMB data are capable to 
put an interesting upper bound to the dark matter - dark energy coupling: 

/3 < 0.16 

(95% c.l.) regardless of the potential (within the class we considered). This implies that the scalar gravity is at least 
1//3 2 w 40 times weaker than ordinary tensor gravity. As shown in Amendola (1999), the limit on (3 can be restated 
as a limit on the constant £ of the non-minimally coupled gravity, £ < 0.017. An experiment like the Planck mission 
can lower the upper bound to /3 to 0.05: scalar gravity would be in this case at least 400 times weaker than ordinary 
tensor gravity. This limit is comparable to those that local gravity experiments impose on the scalar gravity coupling 
to baryons, 0t aryons < 10~ 3 (see e.g. Groom et al. 2000). 

In contrast, CMB data, on their own, cannot put any firm limit to the dark energy potential, unless a narrow prior 
on h is adopted: e.g. h — 0.65 ± 0.05 gives = 0.55 ± 0.2 while h — 0.75 ± 0.05 gives = 0.35 ± 0.2 (assuming 
gaussian priors and marginalizing as usual over all the other parameters, including the coupling). 

In a subsequent work we will consider also the trajectories that do not undergo tracking: in this case the field 
4> is frozen until recently, and the trajectories are indistinguishable from a pure cosmological constant, except for 
the existence of a 0MDE if j3 ^ 0. As expected, we will show that even for off-tracking trajectories the CMB data 
constrain /3 but not the potential. This shows that the constraint on (3 is independent both of the potential and of 
the initial conditions. 
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